/* Ce fichier permet d'évaluer la limite d'un quotient*/
kill(all)$
/* Saisir ci-dessous le numérateur et le dénominateur du quotient */
numerateur:x^2-3*x+2$
denominateur:-3*x^2+5*x+2$
/* Saisir en quel nombre sera effectuée la limite */
limite:2$
/* Saisir 1 pour la limite à droite, -1 pour la limite à gauche */
sens:-1$
/* Saisir le nombre de décimale utilisé dans l'affichage des valeurs */
fpprec:17$
/* Saisir l'emplacement du fichier récupérant les données */
fichier:"/home/user/Bureau/calculInfinetesimal.csv"$
/* ============================================== */
/* Corps du programme */
file_output_append:false$
with_stdout(fichier,printf(true,""))$
file_output_append:true$
f(x):=printf(false,"~h",bfloat(floor(float(x*10^fpprec))/10^fpprec))$
f(5);
a:[]$b:[]$c:[]$d:[]$
for aa:1 thru 8 do 
(xx:limite+sens*10^(-aa),
aaa:float(f(xx)),
a:append(a,[aaa]),
bbb:f(ev(numerateur,x=xx)),
b:append(b,[bbb]),
ccc:f(ev(denominateur,x=xx)),
c:append(c,[ccc]),
ddd:f(ev(numerateur/denominateur,x=xx)),
d:append(d,[ddd]),
print("Pour x=",f(xx),", ",numerateur/denominateur,"=",bbb/ccc,"=",ddd))$
a:simplode(a,";")$
b:simplode(b,";")$
c:simplode(c,";")$
d:simplode(d,";")$
with_stdout(fichier,printf(true,a))$
with_stdout(fichier,printf(true,"~%"))$
with_stdout(fichier,printf(true,b))$
with_stdout(fichier,printf(true,"~%"))$
with_stdout(fichier,printf(true,c))$
with_stdout(fichier,printf(true,"~%"))$
with_stdout(fichier,printf(true,d))$